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Abstract 

The system of two coupled Bose-Einstein condensates is mapped onto a uni- 
axial spin with an applied magnetic field. The mean-field interaction, the 
coupling and the asymmetry or the detuning correspond to the anisotropy, 
the transverse field, and the longitudinal field, respectively. A generalized 
Bloch equation is derived. In the low barrier limit for the quasi-spin model, 
the tunneling rate is analyzed with an imaginary-time path-integral method. 
The dependence of the tunneling rate on the system parameters is obtained. 



* Corresponding author: chlee@mpipks-dresden.mpg.de and chleecn@hotmail.com. 



The crossover temperature T c from the thermal regime to the quantum regime 
is estimated. Below T c quantum tunnelling prevails, otherwise thermal acti- 
vation dominates. 
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I. INTRODUCTION 



The experimental realization of measuring the relative phase and the population oscilla- 
tion between coupled Bose-Einstein condensates (BECs) stimulates great interest in investi- 
gating their macroscopic quantum tunnelling dynamics [1-3]. There are two different types 
of atomic tunnelling between coupled BECs, external tunnelling and internal tunnelling [2,3]. 
The former has different spatially separated single-particle states in a double-well or multi- 
well potential and the latter has different hyperfme internal states in a single-well potential. 
For external tunnelling, the phase interference between BECs confined in a multi-well po- 
tential has been observed [4,5]; the experimental observation of the tunnelling among BECs 
confined in multi-well potential has also been reported [6-8]. For internal tunnelling, J1LA 
realized a two-component BEC coupled with Raman pulses [9] , MIT observed the tunnelling 
across spin domains in BECs [10,11], and LENS reported the current-phase dynamics in two 
weakly coupled BECs trapped in different Zeeman states [12]. 

With the proceeding of the experimental exploration, a lot of theoretical investigation 
was performed simultaneously. Williams et al. demonstrated the existence of Josephson 
tunnelling in a driven two-state single-particle BEC in a single-well trap potential [13]. 
Kasamatsu et al. investigated theoretically the existence of a metastable state and the 
possibility of decay to the ground state through macroscopic quantum tunnelling in two- 
component BECs with repulsive interactions [14]. Smerzi et al. studied the coherent atomic 
tunneling and population oscillations between two zero-temperature BECs confined in a 
double-well potential [15-18]. Macroscopic quantum self-trapping (MQST), namely a self- 
maintained population imbalance with nonzero average value of the fractional population 
imbalance, and 7r— phase oscillations in which the time averaged value of the phase difference 
is equal to ir were detailed in Refs. [15,16]. The authors of Ref. [17] claim that interaction 
with a thermal cloud will damp all different oscillations to the zero-phase mode. In addition, 
macroscopic quantum fluctuations have also been discussed by using second-quantization 
approaches [18,19]. Within the time- dependent potential, chaotic population tunnelling 
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emerges. Abdullaev and Kraenkel analyzed the nonlinear resonances and chaotic oscillations 
of the fractional population imbalance between two coupled BEC's in a double- well trap with 
a time-dependent tunneling amplitude for different damping [20]. They also considered the 
chaotic atomic population resonances and the possibility of stabilization of the unstable- 
mode regime in coupled BEC's with oscillating atomic scattering length [21]. In a previous 
paper, we investigated the chaotic and frequency-locked population oscillation between two 
coupled BECs [22]. 

Although many papers appear in the field of the tunnelling between coupled BECs, 
because of the nonlinearity in the Gross-Pitaevskii equation (GPE), few of them address 
the question of calculating the tunnelling rate and the crossover temperature between dif- 
ferent tunnelling regimes. However, the tunnelling rate and the crossover temperature of 
the spin systems have been studied systematically with the imaginary-time path-integral 
method, including models with applied magnetic field [24-30] and without [31-33]. For a 
two-state system described with linear Schrodinger equation, it is easy to visualize the ef- 
fects of coupling between two states by introducing Bloch's spin vector formalism [23]. Can 
we introduce a generalized Bloch vector for two coupled BECs described with the nonlinear 
Schrodinger equation to map it onto a spin system, and then calculate the tunneling rate and 
the crossover temperature with the imaginary-time path-integral method? If the coupled 
BECs is equivalent to a spin system, the tunnelling process is related to the decay of the 
metastable MQST state to the ground state. More interestingly, the crossover temperature 
corresponds to the transition from the classical or mean-field regime to the second quan- 
tization regime. In the next section, by introducing a generalized Bloch spin vector, the 
coupled BECs is mapped onto an uniaxial spin with an applied magnetic field. In section 
III, the tunneling rate is calculated with the imaginary-time path-integral method, and the 
crossover temperature is estimated. In the last section, a brief discussion and summary is 
given. 
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II. QUASI-SPIN MODEL FOR TWO COUPLED BOSE-EINSTEIN 

CONDENSATES 



Consider the experiments of JILA [9], two Bose-Einstein condensates in the \F = l,m F = 
— 1 >= |1 > and \F = 2,m F = 1 >= |2 > spin states of 87 Rb are coupled by a two-photon 
pulse with the two-photon Rabi- frequency Q and a finite detuning 5 = uja — oJhf- Where, 
ujd = u>i +u 2 is the driven frequency of the two-photon pulses, Uhf is the transition frequency 
between two hyperfine states. In the rotating frame, ignoring the damping and the finite- 
temperature effects, the coupled two-component BEC system can be described by a pair of 
coupled GPEs 

in ^A = {H o + h mf _ |^ 2( - t) + t)f 
ih 9^A = {H o + h mf + [f^^r, t) + f * 2 (r , t), 

where, the free evolution Hamiltonians H® = + Vi(r) (i = 1,2) and the mean- field 

interaction Hamiltonians Ef F = ^(a;;|^(r , t)\ 2 + a^jfj , t)\ 2 ) = 1,2, % ^ j). The 
coefficient is the scattering length between states i and j and it satisfies = dji. Weak 
coupling is defined by the Rabi frequency satisfying Q/ 1 (uOxUOyUOz) 1 ^ = Q/cJ 1, where 
ZJ = (u x ujyUJ z y/ 3 is the geometric-averaged angular frequency for the trapping potential. 
In this regime, we can write the macroscopic wavefunctions using the variational ansatz 
Vi(r,t) = ipi(t)®i(r) with ip^t) = y/Ni(t)e iai ^ (i = 1,2). In the ansatz, the functions 
<3>i(r) describe the spatial distribution of the % — th component, the complex coefficient 
functions ipi(t) are spatially uniform and contain all time-dependence in the macroscopic 
quantum wave- functions ^i(r,t). The symbols Ni(t) and ««(£) represent the populations 
and phases of the i — th condensate, respectively. Because the coupling is very weak, the 
spatial distributions vary slowly in time and are very close to the adiabatic solutions to 
the time-independent uncoupled case for GP equations (1), being slaved by the populations 
[13]. Thus, the complex coefficient functions ipi(t) obey the nonlinear two-mode dynamical 
equations 
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ihi^it) = (E° 2 - f + u 22 \^(t)\ 2 + u 21 \Mt)\ 2 )Mt) + fA(t), 

ihf^it) = (£° + f + C/nl^WI 2 + tfi 2 |& (*)!%(*) + f 

The parameters satisfy £° = / r r)d r, C^- = / |^(r )| 2 |$ i (r )| 2 d r= [7^ 

and K = hVt j $>i(r)Q 2 (r)d r = 1,2). The terms in K describe population transfer 
(internal tunnelling) between two BEC states, whereas the terms in Uij, which depend 
on the numbers of atoms in each BEC state, describe the mean-field interaction between 
atoms. When U 2 i and 5 equal zero, these coupled equations can also describe the BECs 
in a double- well potential [15-18]. Similar to the coupled two-state system obeying the 
linear Schrodinger equation, we introduce a generalized Bloch spin vector (u, v, w) with the 
components 

u = Mi + Ml » v = -KMl - Mi) > w = M* - Mi- ( 3 ) 

Obviously, u 2 + v 2 + w 2 = (Aq + N 2 ) 2 = N% is a conserved quantity when finite-temperature 
and damping effects can be ignored. Rescaling the time t/h to t, the Bloch spin vector 
satisfies 

du , . dv T ^ . . dw , A , 

— = u(7 + r]w) , — = Kw- u(7 + rjw) , — = -Kv, (4) 

where 7 = £ 2 ° - £7° + iV T (f/ 22 - C/n)/2 - HS and 77 = (U 22 + E/n - 2U 12 )/2. Regarding 
the atom in one condensate as spin-up state and the atom in the other condensate as spin- 
down state, the coupled BECs can be described with the quasi-spin S— u e x +v e y +w e z . 
In this language, the longitudinal component w depicts the population difference, and the 
transverse components u and v characterize the coherence. Thus the effective Hamiltonian 
for the quasi-spin is 

E = - 1 - V S 2 Z -KS X - 1 S Z . (5) 

The above Hamiltonian is similar to the one of a uniaxial spin with an applied magnetic field 
[26-30], it indicates that the mean-field interaction brings the anisotropy 77, the coupling 
causes an effective transverse magnetic field K along axis-x, and the asymmetry or the 
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detuning induces an effective longitudinal magnetic field 7. In the symmetric case (E® = E®, 
U22 = Uu and 8 — 0), it is consistent with the one derived from the second quantized 
Hamiltonian in [4]. 

III. TUNNELLING RATE AND CROSSOVER TEMPERATURE 

In conventional spherical coordinates, the spin components can be written as S x = 
N T sin cos (f), S y = N T sin 9 sin <p and S z = N T cos9 (see Fig. 1). Thus, the corresponding 
effective Hamiltonian is formulated as 

IK 

E = -riNri- cos 2 9 + sin# cos0 + —— cos 9). (6) 

' TV 2 r]N T 7]N T 

Based upon the analysis of a spin in a uniaxial magnetic field [26-30], we know that there 
are stationary states if some angles (# o ,0o) satisfy dE/d<p\ e ^^ Q = and dE/d9\ r p:£ o = 0. 
The condition dE / <90|^, o o = locates the stationary states in the XOZ plane (sin0 o = 0). 
The existence of multiple stationary states in this quasi-spin system is equivalent to the 
existence of multiple metastable MQST states in the coupled BECs. Near the metastable 
states the potential describes a "canyon" satisfying 

E e = E(9, O )/(^|) = ~\ c o s2 e ~ p co *( d ~ Op)- (7) 

The parameters obey P = \JK 2 + ^y 2 /\r]N T \, sin9p = K cos O / V K 2 + 7 2 and cos9p = 
7/a/-^ 2 + 7 2 . As stated in the previous section, the parameter K oc Q > 0, therefore 
sin#p > and sin#p < correspond to the equal-phase mode (0o = 0) and the anti- 
phase mode (0o = ^) in the coupled two-component BECs, respectively. In the case of 
E® — E\ + N T (U 2 2 — Un)/2 = 0, the parameter 7 is just the negative detuning —5, thus 
cos9p > and cos9 P < correspond to the red detuning and the blue detuning of the 
coupling laser, respectively. The dE /d9\^l^ o = is equivalent to dEe/ d9\ e ^^ Q = 0, that is, 
sin 2^o + 2Psin(# — Op) = 0. For some critical points where both the first and the second 
derivatives of Eg equal zero, an appreciable tunnelling rate appears. This gives 
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sin 29 c + 2P C sm(9 c - 9 P ) = 0, 

( 8 ) 

cos 29 c + P c cos(9 c - 9 P ) = 0, 
here, 9c and Pc are critical values for 9 and P, respectively. Solving the above equations, 
one can obtain tan 3 9 C = — tan6> P and P c = (sm 2 ^ 3 9 P + cos 2 / 3 9 P )~ 3 / 2 . The system has an 
instanton solution at the critical point P = P c , i.e., (K/(r]N T )) 2 / 3 + ( , y/(r]N T ) ) 2 / 3 = 1. This 
critical point stands on the separatrix between the single stable regime and the multiple- 
stable regime. It separates the metastable multi-MQST behavior between the single-stable 
population oscillation in the coupled two-component BEC. 

According to the dependence of E e on 9, we obtain that the condition for the existence of 
multiple stationary states is P < P c , i.e., (K/(r)N T )) 2 / 3 + (7 /{i]N T )) 2 / 3 < 1. One can easily 
find the small oscillations around these stationary states with nonzero time-averaged values 
for S z and ^S 2 + S 2 . These oscillations correspond to the phase-locked MQST states with 
time-averaged relative phase or n and multiple stationary states correspond to multiple 
metastable MQST states with fixed nonzero population difference and relative phase or 
7r. The appearance of multiple stationary sates indicates, only for some proper parameters, 
that multiple metastable MQST states exist. For simplicity we only consider the case where 
the parameter P is slightly lower than the fixed critical value Pc, P = Pc(l — e), e 1. 
This requires that the Rabi frequency, the detuning, the scattering lengths, and the total 
atomic number in the coupled BEC system must cooperate with each other to approach the 
critical values for the emergence of multiple metastable MQST states. One way to maintain 
the critical value Pc unchanged is fixing the values of the ratio 7/ if and other correlated 
parameters (77 and N p ), that is, keeping the angle 9 P unchanged. By introducing a new 
positive variable £ = 9 — 9 , the potential (6) can be expanded into 

E e {9) = E e (9 ) + i[Vfe£ 2 - £ 3 + 0(£ 4 )] sin(2£ c ). (9) 

With the definition in Refs. [26,27,29,30], the tunnelling rate T obeys N p (t) = 
iVp(O) exp(— Ft) and it can be written as T = Aexp(-B) for the quantum tunnelling 
regime. Here, N P (t) is the population occupying the metastable state at time t and the 
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tunnelling exponent B (> 0) is determined by the imaginary time action of the instan- 
ton solution. Similar to Ref. [26], the tunnelling exponent follows from the path integral 
f _D{0(r)} f D{cos 0(t)} exp(I/h) over the continuum of trajectories which start and end 
at (#o? 0o ) an d which are close to the instanton solution, where, r is the imaginary time it, 
and / is the imaginary time action I — J dr[iN T (l — cos 9) dtp /dr + E(0, 0)]. Integrating the 
imaginary time action by parts, one can gain the tunnelling exponent 



B = 16 x 6 1/4 iV T [l - (1 + | 7 /^D(1 + \l/K\ 2 / z )-^\K/{r ] N T )\f A \ 1 /K\ 1 / & /b. (12) 



To control the tunnelling, one has to select proper values for parameters 7, K and rj. In 
the experiments performed in a double- well potential [4,5], it can be realized by modifying 
the barrier position, the barrier height and the magnetic field (using Feshbach resonances 
to adjust the scattering lengths, [35]), respectively. In the experiments with two-component 
BECs in a single-well potential [9] , it can be realized by adjusting the laser detuning, the 
laser intensity, and the magnetic field, respectively. For fixed value of rj and j/K, the 
tunnelling exponent B decreases with the increasing of the intensity of the coupling laser. 
In Fig. 2, we show how the tunnelling exponent B depends on the angle Op. In the region 
between and 7r, the ratio B(0 P ) / 'B(n /4) decreases from positive infinity to zero when the 
angle Op equals n/2, which corresponds to the symmetric case (7 = 0), and then increases 
to positive infinity when the angle Op is close to it. It is almost flat when the angle Op is 
not close to 0, rr/2 and n. This angular dependence indicates, in the case of fixed value of 
e, that the tunnelling exponent increases with increasing \ j/K\. 




(10) 



From the definition of e, one can obtain 



e = l- P/P c = 1 - (1 + |7/^| 2 )(1 + ll/Kl^r^lK/^Nr)]. 



(11) 



Thus the tunnelling exponent can be expressed as 
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The result for the angle 6 P close to n/2, which corresponds to the symmetric case 7 = 0, 
should be taken with great caution because the coefficient sin(26> c ) in the Taylor expansion 
series (9) is equal to zero. In this case, the problem corresponds to the tunnelling between 
two equivalent minima which correspond to the angle 6p equal and n. Thus the potential 
can be expanded into the form of £ 2 — £ 4 and the tunnelling exponent B is expressed as 
B = ASs 3 / 2 = ANt€ 3 ^ 2 . Therefore, the tunnelling exponents (10) and (12) only hold for the 
asymmetric case where 7 7^ 0. 

To confirm our prediction from the quasi-spin model, we perform a numerical simula- 
tion of the equation (2). A qualitative change in the stationary-state behavior occurs at 
\K/(rjN T )\ = 1. When \K/(rjN T ) | > 1, there are no metastable states for any effective 
detuning 7. However, when \K/(rjN T )\ < 1, metastable states exist in the region [— 7 C , +7J 
for proper relative phase, where 7 C satisfies (K/(r]N T )) 2 ^ 3 + ('yj \r]N T )) 2 / 3 = 1. See the left 
column of Fig. 3. Two stationary states, indicated as S\ and S 2 in the figure, are stable 
and the other one (U) is unstable. Adiabatically changing the effective detuning 7 from 
7 C — e to 7 C + e (e is a very small positive number), in the space of the fractional population 
difference z = (N 2 — Ni)/N T and the relative phase = a 2 — cxi, a trajectory in the vicin- 
ity of £2 becomes a large orbit C encircling Si. From the views of instanton method, the 
tunnelling exponent is determined by the canonical action of the orbit, i.e., B follows from 
the path integral J D{z(r)} J D{<p{r)} exp(I c /h) over the continuum of trajectories which 
are close to the instanton solution. At different bifurcation points j c , the numerical results 
show B(\-f/K\)/B(\-f/K\ = 1) oc | 7 /jqo.i63±o.oo2 w ^/k] 1 / 6 , this confirms our previous 
prediction from the quasi-spin model (see the right column of Fig. 3). 

There are two important aspects which must be noted. The one is that these results for 
tunnelling are only valid in the low barrier limit for the quasi-spin model, i.e., e « 1. This 
means that the above results only hold in the region which approaches the critical point of 
emergence of multiple metastable MQST states. The parametric dependence of the general 
case is still an open problem. The other is the validity of the Wentzel-Kramers-Brillouin 
(WKB) semiclassical approximation. The semiclassical approach can only be used in the 
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case of small tunnelling probability, that is, B » 1. In this low barrier limit, from the 
Taylor expansion series (9) one can obtain the following tunnelling amplitude by using the 
theory developed by Caldeira and Leggett [34], 

A = {IbB/^f^uj, 

= ^ T (l|5)l/2(|)l/4| CQt p |l/6/(! + cot 2/3 Q p ^ (13) 
= ^ T (^)V2(|)V4| 7/K |l/6 /(l + | 7 /K|2/3). 

Here, u is the angular frequency of small oscillations near the bottom of the inverse potential. 
Apparently, when the angle Op is close to the kir (k = 0,1), which corresponds to small 
Rabi frequency or large detuning of the coupling laser between two BECs, the tunnelling 
amplitude A approaches to zero, see Fig. 4. As presented in above, in the case of Op close to 
7r/2 which corresponds to the symmetric case 7 = 0, the potential is not in form of £ 2 — £ 3 
but in form of £ 2 — £ 4 because the coefficient $m(2Q c ) in the Taylor expansion series (9) 
equals zero. Therefore, the above formula for the tunnelling amplitude only holds for the 
asymmetric case 7^0. Generally, contrary to the tunnelling exponent B, the tunnelling 
amplitude A is sensitive to the structure of quantum levels in the potential. Therefore, for 
the case of the full potential (6) and (7), the estimation of A is still an open problem. 

Population transfer between two states in a bistable system can occur either due to 
classical thermal activation which depends on the system temperature or due to quantum 
tunnelling which does not depend on the system temperature. There exists a phase transition 
from the thermal regime to the quantum regime which occurs at the crossover temperature 
T c . Above Tc, quantum effects are very small and the population transfer rate follows the 
Arrhenius law, 

U B 

r^ erma « = r exp(-^-^). (14) 

Here, Ub is the height of the energy barrier between two states and ks is the Boltzmann 
constant. Below Tc, the population transfer is purely quantum, 

T quantum = Aexp(-B), (15) 
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with B independent of the system temperature. Thus the transition occurs when Thermal = 
^quantum- Neglecting the prefactors and equating the exponents, the crossover temperature 
can be estimated as 

T c = U B /(k B B). (16) 

The transition region is approximately the temperature interval [T c (l — B^ 1 ), T c (l + B~ L )]. 
This crossover resembles a first-order phase transition of the tunnelling rate T because it is 
accompanied with the discontinuity of dT/dT at T c [29]. 

There is another regime for tunnelling, the thermally assisted tunnelling (TAT), in which 
the particle strides over the barrier to the bottom of the potential with lowering temperature 
[29,30]. The transition from the classical regime to the TAT regime resembles a second-order 
classical-quantum phase transition of the tunnelling rate T because it is accompanied with a 
discontinuity of d 2 T/dT 2 and no discontinuity of dT/dT at the crossover temperature. The 
corresponding transition temperature can be estimated as 

T> c = h/(r k B ) = fuj/(2nk B ), (17) 

where r and oj are the period and the angular frequency of small oscillations near the bottom 
of the inverse potential, respectively [25-27,29,30]. In the low barrier limit (e « 1), from 
the Taylor expansion series (9), one can obtain the barrier height 

U B = 77iV|(2 £ /3) 3/2 | sin(20 c )| = ^(h/r )B, (18) 

where, 

|sin(2^)|=2| 7 /^ c /(l + | 7 /i^). (19) 

Comparing both crossover temperatures, one can easily find that they differ by a factor 
T c /T[, = 5tt/18 = 1/1.15, which means that they are of the same order of magnitude and 
can both be used to estimate the crossover temperature. 

Below, from the experimental parameters in the experiments of JILA [9], we will give 
a quantitative estimation for the tunnelling rate and the crossover temperature. In those 
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experiments, the atomic mass mi = m 2 = = 1.45 x 10 25 kg, the time- averaged orbiting 
potential (TOP) magnetic trap has an axial frequency v z = 59 Hz and a radial frequency 



v z /V& = 21 Hz, the s-wave scatter lengths an = 5.36 nm, au = &2i 



5.53 nm and a 2 2 = 5.70 nm, and the total atomic number N T «5x 10 5 . To obtain the 
numerical values conveniently, we choose the natural units of the problem, in which, time 
is in units of l/(u x u y u} z y/ 3 = 1/cJ, length is in units of the size of the geometric-averaged 
harmonic-oscillator length d = -\/h/[ {uj x uj y uj z ) ^mm] = a/ h/ (Umjn,), energy is in units of 
the geometric-averaged trap level spacing h{uo x uo y uo z ) 1 ^ = KlJ, and mass is in units of Rb 
atomic mass mja- 

Due to gravity acting besides the TOP, the centers of two condensates will displace 
along the vertical direction and the two equilibrium displacements are generally not the 
same. Thus, if the interparticle interaction is absent, the lowest single-particle state has the 
familiar wave function, 

*»-< ? » - exp( -4 - $ - ^w^- (2o » 

Where, F i = +1, F2 = — 1, 2^ is the offset between two potential centers along the vertical 
axis, dk = y/h/ (uJk m Rb) (k = x, y, z) are the oscillator lengths. The offset 2z between two 
condensates can be varied by adjusting the magnitude of the rotating magnetic field. In 
the presence of interatomic interaction, the dimensions of the condensates are changed. The 
spatial parts of the macroscopic quantum wave functions are in the shape of 

*- (? > = ^w„y» ^'i - W, ~ {j ^t ly (21) 

The variational parameters bit (k — x,y, z; i — 1,2) depend on the scattering length, the 
total atom number, and the trapping potential and they have almost the same numerical 
values as aV For proper values of the offset 2z , the numerical results of [13] show that the 
spatial distributions $j(r) and their overlap only weakly depend on the total atom numbers 
in each condensate. For simplicity, in the following calculations, the variational parameters 
bit are replaced by the oscillator lengths dk- Therefore, the parameters Ef, Uij and K are 
determined by 
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El = El = h(u x + u y + u z )/2, 

Uu = 4n h 2 aii/[(V^) 3 d x d y d z m Rb },(i = 1,2), 

(22) 

U12 = 4nh 2 a 12 exp(-2z%/dl)/[(V2lr) 3 d x d y d z m Rb ], 
K =hnexp(-z 2 /d 2 z ). 

So the corresponding parameters in the quasi-spin model (5) can be written 
as 7 = h 2 N T (a 22 — flu) / \V2rrd x d y d z mRb) — h5 and rj = h 2 [a 22 + oil — 
2a\2 exp(— 2z 2 /d 2 z )}/(V27Td x d y d z m Rb ). In the case of complete overlap (2z = 0), the 
anisotropy parameter rj equals zero, thus the metastable multi-MQST behavior will never 
appear, but some running-phase MQST states may still exist. This indicates that, to 
insure the existence of multiple metastable MQST states, a finite offset must be kept 
between two condensates. Furthermore, the appearance of this kind of MQST requires 
K 2/s + 7 2/s < ( vNT y/\ Because K oc Q and 7 oc 5, this inequality indicates that the Rabi 
frequency and the detuning of the coupling pulses must be relatively small. Choosing the to- 
tal atom number iV T = 2.0 x 10 4 , the half offset z = 0.20d 2 , the Rabi frequency f2 = 2ir x 10 
Hz, and the detuning 5 = -179 Hz, one can get 7]N T = 6.70 x 10 -32 , 7 = 4.57 x 10~ 32 , 
K = 6.90 x 10~ 33 and e = 9.78 x 10~ 3 . Thus, the corresponding tunnelling exponent B and 
crossover temperature T c are around 4.22 x 10 2 and 3.54 x 10~ 2 nK, respectively. Obvi- 
ously, the crossover temperature Tc, which corresponds to a phase transition from classical 
tunnelling to quantum tunnelling, is far below than the critical temperature T ~ 150 nK 
for Bose-Einstein condensation in a dilute gas of 87 Rb. 



IV. DISCUSSION AND SUMMARY 

The generalized Bloch equation (4) and its stability analysis will help to control the pop- 
ulation transfer and realize the single-qubit operation with BECs qubit. Theoretically, any 
two-state quantum system can serve as a qubit, many of them have been realized experimen- 
tally. To make use of two quantum states, the coherence and superposition between them is 
the most essential qualification. The experimental observation of coherence and superposi- 
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tion between two BECs indicates the possibility of encoding two coupled BECs as a qubit. 
However, because of the mean-field interaction among Bosonic condensed atoms, the qubit 
operations become very difficult to perform. To accomplish a single-qubit operation, it must 
be possible rotated arbitrarily in the Hilbert space. This requires the atomic populations can 
be transferred arbitrarily. From the Bloch equations (4), we find that MQST prevents the 
arbitrary rotation of the state vector. And even if there no MQST, when r\ ^ 0, the complete 
population inversion can not be accomplished with linear operations. Thus, to accomplish a 
linear qubit operation, one has to adjust the parameter r\ to zero by varying the atomic scat- 
tering length with a Feshbach resonance [35]. In this case, the mean-field interaction gives 
a density-shift to the original energy levels and, according to Rabi's theory, the arbitrary 
rotation of the state vector can be performed easily. Thus, if one encodes the qubit states 
|0) and |1) as the condensate wavefunctions for two condensates in a double-well potential or 
two hyperfine-state condensates coupled with Raman pulses [36] , an arbitrary one-bit linear 
operation can be realized when the anisotropy is absent [r\ = 0) and an arbitrary one-bit 
nonlinear operation can be realized when the metastable multi-MQST behavior is absent 
(\K\ > I^A^I). This means that, to perform an arbitrary one-bit transformation, it at least 
needs choosing proper parameters to avoid the emergence of the metastable multi-MQST 
behavior. 

The tunnelling of the quasi-spin model described by the Hamiltonian (5) has also been 
investigated by mapping it onto a particle moving in an asymmetric double-well potential 
[27-29]. Using this approach, Garanin et. al. have explored some new fascinating feature 
of this uniaxial spin model in the strongly biased limit [29]. They find that there exist 
two different regimes for the classical-quantum transition of the tunnelling rate and the 
kind of transition depends on both the strength and the direction of the magnetic field. 
In this article, we directly analyze the tunnelling in the low barrier limit for the quasi- 
spin model, which corresponds to the effective magnetic fields near their critical values for 
appearance of metastable states. This requires that all physical parameters of the coupled 
BECs collaborate with each other to approach the critical point of appearance of multiple 
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metastable MQST states. The symmetric case (7 = 0) of the coupled BECs corresponds to 
the unbiased case (H z = 0) of the anisotropic spin model, which has been investigated in 
details by mapping it onto a particle moving in a symmetric double-well potential [25]. 

The macroscopic quantum tunnelling of two-component BECs has also been investigated 
by the Kasamatsu group. Using a numerical approach, they have analyzed the tunnelling 
between two kinds of metastable stationary states, a symmetry-breaking state (SBS) and a 
symmetry-preserving state (SPS), in uncoupled two-component BECs [14]. To improve the 
usual Gaussian variational method, they have introduced a collective coordinate approach 
and then calculated the tunnelling rate within the WKB approximation. In that system, 
the populations of the two components can not be converted into each other because of the 
absence of coupling. This means, the tunneling does not occur between two components but 
between stationary states with different spatial configurations. Thus, this kind of tunnelling 
originates from the quantized spatial structure of the Hamiltonian. In our model, due to the 
coupling, the population can be transferred from one component to the other. Furthermore, 
we assume the coupling is very weak, thus both components stay in their ground stationary 
states through the full process. The metastability (metastable MQST) is the result of the 
cooperation between the coupling and the mean-field interaction (including both the intra- 
component and the inter-component interaction). Correspondingly, the tunneling from the 
metastable self-trapped state to its ground state of the coupled two-component BECs is 
caused by the quantized structure of their field operators. 

In conclusion a system of coupled BECs (two BECs in a double-well potential or two 
internal state BECs coupled with laser pulses) has been mapped to a spin in a magnetic field 
by introducing a generalized Bloch vector. The mean-field interaction, the coupling and the 
asymmetry or the detuning are relevant to the anisotropy, the transverse magnetic field and 
the longitudinal magnetic field, respectively. The corresponding generalized Bloch equation 
is obtained. The analysis of this generalized Bloch equation will be propitious to control the 
population transfer and realize the quantum computation with coupled BECs. Based upon 
experience from the well-studied tunneling of spin systems, the detailed information about 
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the tunnelling between two metastable MQST states in coupled two-component BECs can be 
obtained with the imaginary-time path-integral method. The crossover temperature T c at 
the critical point for a transition from the classical thermal regime to the quantum regime was 
obtained. When the system temperature decreases through Tc, the population conversion 
goes from classical thermal activation regime to purely quantum tunnelling regime. This 
means, below the crossover temperature Tc, the quantum fluctuations in the atomic fields 
take the dominant position. We also find that the tunnelling rate can be adjusted by varying 
the coupling and the trapping magnetic field. 
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Figure caption 

Fig. 1 The quai-spin S and it's components (u, v, w) in conventional spherical 
coordinates. 

Fig. 2 The tunnelling exponent ratio B{0p) / B{0p = 7r/4) versus different 
9p. Where, the angle P characterizes the angle between the effective magnetic 
field B e ff= K e x +7 e z and axis-z. 

Fig. 3 In the left column, the stationary states for \K\ < \i]Np\ are shown. 
There are two metastable states Si, S 2 and one unstable state U. In the right 
column, the tunnelling exponent ratio B/B versus different \j/K\ is presented, 
where B = B{\^/K\ = 1). The black dots show the numerical data and the 
straight line represents the linear fit for the logarithmic data. 

Fig. 4 The tunnelling amplitude ratio A(0p)/A(0 P = 7r/4) versus different 
Op, where the angle Op characterizes the angle between the effective magnetic 
field B e ff= K e x +7 e z and axis-z. 
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